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Abstract 

We study a two-dimensional fluid of dipolar hard disks by Monte Carlo simulations in a square 
with periodic boundary conditions and on the surface of a sphere. The theory of the dielectric 
constant and the asymptotic behaviour of the equilibrium pair correlation function in the fluid 
phase is derived for both geometries. After having established the equivalence of the two methods 
we study the stability of the liquid phase in the canonical ensemble. We give evidence of a phase 
made of living polymers at low temperatures and provide a tentative phase diagram. 
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I. INTRODUCTION 


This paper is devoted to a study of a two-dimensional (2D) system made of identical 
dipolar hard disks (DHD) in the Euclidian plane E 2 by means of Monte-Carlo (MC) sim¬ 
ulations. The dipoles are assumed to be permanent and the configurational energy of N 
dipolar molecules in E 2 reads as 
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In Eq. ([T]), Vusi^) is the hard disk potential of diameter a. The second term is the contri¬ 
bution from the 2D dipole-dipole interaction where /Xj = /is,, /i permanent dipole moment, 
Sj unit vector in the direction of the dipole moment of particle i, = rj — Uj, the vector 
joining the centres of mass of the particles, and Vij = We stress that the system that 
we consider cannot be seen as a thin layer of a real 3D system of dipoles. In this case, the 
electrostatic interactions should be derived from the solutions of the 3D Laplace equation 
while the dipole-dipole interaction involved in Eq. ([T]) is derived from the solution of the 2D 
Laplace equation in the plane. 

We have performed MC simulations of the DHD fluid in a square with periodic bound¬ 
ary conditions and on the surface of an ordinary sphere. In both cases the dipole-dipole 
interaction is obtained from a rigorous solution of Laplace equation in the considered ge¬ 
ometry [l-4|. We compare the two methods in the liquid phase and check that they both 
yield the same thermodynamic, structural and dielectric properties. Both methods are then 
used for preliminaries MC studies of the DHD fluid at low temperatures. In this regime, 
as for real 3D dipoles confined in a plane (see e.g., Ref. j^), the 2D dipoles aggregate to 
form living chains and ring polymers at low densities and more involved structures at higher 
densities. 

The paper is organized as follows. After this introduction we give details on the two 
simulation techniques used in this work in Sec. Ull Next Sec. UlI] is devoted to a digest of 
the general theory of dielectric media in an arbitrary 2D geometry jj], with applications to 
the square with periodic boundary conditions and the sphere. This theoretical analysis is 
notably required to understand the long range tails of the pair correlation functions in both 
geometries. Checks of these asymptotic behaviours as well as quantitative comparisons 


2 








between the two methods are discussed in Sec. m In Sec. m we present extensive MC 
simulations of the DHD fluid by both methods and give a tentative phase diagram of the 
system. We conclude in Sec. El] 


II. SIMULATION METHODS 


A. Periodic boundary conditions 


1 is a square of side L with periodic boundary conditions, 
, ^ . Some care is required to take into account the long 


In this method the simulation ce 
that will be referred to as space C 2 
range of dipole-dipole interaction. The usual way to compute the conhgurational energy Um 
is to replicate the basic simulation cell periodically in space and calculate Udd as the sum of 
the interactions of the N dipoles in the basic cell with all the other dipoles in the cell and 
with the periodically repeated images in the surrounding cells 
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The prime affixed to the sum over n = {nx,ny), with rix^Uy integers, means that the term 
i ^ j is omitted when n = 0. 

By a lattice summation technique (Ewald sum) the slowly and conditionally convergent 
sum is transformed into two rapidly convergent sums, one in direct space, the other in 
reciprocal space, the rate of convergence of both sums being regulated by the parameter a. 
The resulting expression for the energy of the 2D system is [l| 
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where the functions b{r) and c(r) are given by 
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and 


N 

F{k) = ^ (k ■ Si) exp[ik ■ r*]. (6) 

i=l 

In Eq. ([3]) S' = is the area of the simnlation cell, N the nnniber of particles and F* the 
complex conjngate of F. The wave-vectors k which enter the reciprocal space contribntions 
to the energy are of the form 

k = 27rn/L . (7) 

Care has to be taken to properly choose the a parameter which governs the rate of 
convergence of the real- and reciprocal-space contribntions in Eq. ([3]). It is generally taken 
snfliciently large so that only the terms with n = 0 need to be retained in Eqs. ([3]). The 
last term in Eq. (|3]) represents the contribntion to the energy from the depolarization held 
created when a continnons medinm of dielectric constant e' surronnds a disk shaped sample 
of periodic replica. For a condncting medinm (e' = oo) this term vanishes while for a system 
in vacnnm (e' = 1) it is where M is the total polarization of the system. 

A thermodynamic state of the DHD hnid is characterized by a rednced density p* = 
Na'^/S where S = is the snrface of the sqnare of simnlation and the rednced dipole p* 
with /i*^ = /{ksTa'^) {ks Boltzmann constant, T temperatnre). 


B. Spherical boundary conditions 


In this method the simnlation cell is the snrface of an ordinary sphere of center O and 
radins R, that will be referred to as space S 2 [31,14. The electrostatics can be solved exactly 
in 1 S 2 in two different ways and therefore two distinct models are available { 4 , 

In the first version the DHD flnid is made of N ordinary (or mono-) dipoles = pSi 
tangent to the sphere S 2 at points OM* = Rzi (z* ■ Sj = 0). In the second version considered 
in this article, one rather considers a collection of N bi-dipoles. A bi-dipole is dehned as 
a dnmbell of two identical mono-dipoles located at two antipodal points of the sphere at 
points OMj = Rzi and OMj = —i?Zj. The nnmerical experiments of Ref. jj] show that the 
convergence to the thermodynamic limit is in general faster for bi-dipoles than for mono¬ 
dipoles. The conhgnrational energy of the DHD flnid reads 
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where r’Hs('^*i) hard-core pair potential dehned by 
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where V’p is the angle between vectors Zj and Zj, i.e. cos'i/jij = Zj, -Zj and thus Vij = R^|Jij is 
the length of the geodesic length between points Mi and Mj. The dipole-dipole interaction 
is given by 
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In Eq. dH]) the vectors Zj can always be chosen in the northern hemisphere S 2 because of 
the special symmetries of the interaction. It is thus clear that the actual domain occupied 
by the fluid is the northern hemisphere S 2 rather than the whole hypersphere. In terms of 
mono-dipoles the interpretation of the model is therefore the following : when a mono-dipole 
leaves S 2 at some point Mj of the equator the same dipole moment /Xj reenters S 2 at the 
antipodal point Mj. Therefore bi-dipoles living on the whole sphere are equivalent to mono¬ 
dipoles living on the northern hemisphere but with special boundary conditions ensuring 
homogeneity and isotropy at equilibrium (in the case of a fluid phase). We stress that the 
expression ffTOj) has been deduced rigorously from the solution of Laplace-Betrami equation 
in 1 S 2 jd] by contrast with the heuristic dipole-dipole interaction used in reference j^. 

A thermodynamic state of this model is now characterized by a dimensionless number 
density p* = Na'^/S where S = 271R^ is the 2D surface of the northern hemisphere S 2 and 
the reduced dipole p* with = p'^KhsTa'^). 


III. FULTON’S THEORY 

Let us consider quite generally a polar fluid occupying a 2D surface A with boundaries 
dA. We assume the system to be at thermal equilibrium in a homogeneous and isotropic 
fluid phase. The fluid behaves macroscopically as a dielectric medium characterized by a 
scalar dielectric constant e. Due to the lack of screening in such fluids, the asymptotic 
behaviour of the pair correlation function is long ranged and depends on the geometry of 
the system, i.e. its shape, size, and the properties imposed to the electric held (or potential) 
on the boundaries dA as well. As a consequence, the expression of the dielectric constant 
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e in terms of the fluctuations of polarization also depends on the geome 


ry. These issues 


can be formally taken into account in the framework of Fulton’s theory which achieves 
an elegant synthesis between the linear response theory and the electrostatics of continuous 
media. 

In addition to provide an expression for the dielectric constant e Fulton’s formalism also 
yields the asymptotic behaviour of the pair correlation function. Fulton’s formalism can be 
extended without more ado to non-euclidian geometries and was applied notably to 3D cubic 
systems with periodic boundary conditions (space C 3 ) and hyperspheres 1 S 3 in Refs. BQ, 
and, recently, to the 2D euclidian plane E 2 and the sphere 1 S 2 , for both mono- and bi¬ 
dipoles In this section we derive the missing results for space C 2 and recall the results 
for 2D polar fluids in E 2 and 82 - 

Fulton’s relations constitute the quintessence of Fulton’s formalism; they are formally 
independent of the geometry and read 


X =(J -|- cr o G o cr 

G =Go o (I - 0-0 Go)"^ 


( ll a) 

( ll b) 


Some comments seem appropriate. Let us hrst dehne the tensorial susceptibility X- Under 
the influence of an external electrostatic held S{y) the medium acquires a macroscopic 
polarization 

P(r) =< P(r) >£ , (12) 

where the brackets denote the equilibrium average of the microscopic polarization P(r) = 
Yj) in the presence of the external held S. The relation between the macro¬ 
scopic polarization P and the external held £ can be established in the framework of linear- 
response theory, provided that £ is small enough, with the result 

27rP(ri) = [x o £] (n) j x(ri, ra) ■ ^(rs)^ . (13) 

The r.h.s. of Eq. has been formulated in a compact, albeit convenient notation that 
will be adopted henceforth, where the symbol o (which also enters Eqs. flTT|) l means both 
a tensorial contraction (denoted by the dot ” • ”) and a spacial convolution over the whole 
domain A hlled by the medium. The tensorial susceptibility x Eq. (1 11 ah reads 

X(ri,ra) = 27r/3 < P(ri)P(ra) > , (14) 
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where /3 = l/ksT and the thermal averages are computed in the absence of the external 

held e. 

However, the dielectric properties of the huid are characterized by the dielectric tensor e. 
In Eqs flTT]) we have introduced, following Fulton, the convenient notation cr = e — I with 
I(ri,r 2 ) = U5‘'^^(ri2) where U = e^-ea, + &y&y is the unit dyadic tensor. The tensor cr enters 
the constitutive relation 

27rP = cr o E , (15) 

where the Maxwell held E(r) is the sum of the external held ^(r) and the electric held 
created by the macroscopic polarization of the huid. Therefore one has 

E = 5 + 27rGo o P , (16) 

where Go denotes the bare dipolar Green’s function. Note that 27rGo(ri,r2) ■ ^^2 is the 
electric held at point ri created by a point dipole /X 2 located at point r 2 in vacuum and in 
the presence of the boundary dK. In the presence of the dielectric medium this held is now 
given by 27rG(ri,r2) ■ ^^2 where the macroscopic, or dressed. Green’s function is given by 
Eq. flllb|l in which the inverse must be understood in the sense of operators. 

It is generally assumed that e is a local function, i.e. e = el. More precisely, it is 
plausible -and we shall take it for granted- that e(ri,r 2 ) is a short range function of the 
distance between the two points ri and r 2 , at least for a homogeneous liquid (or in the bulk 
in the presence of interfaces), and one then dehnes 

e\J = [ (fr 2 e(ri,r 2 ) . (17) 

Ja 

Experiments show, and this fact must be admitted, that while e is an intrinsic property of 
the medium Eqs. (ITT]) show that the susceptibility tensor x(ri,r 2 ) is a long range function 
of ri 2 which depends on the considered geometry. The locality assumption on e allows an 
explicit calculation of the Green’s function G(ri, r 2 ) in some geometries, notably those used 
in MG simulations. 

A. The square C 2 

In Ref. Fulton’s formalism was applied to the 3D dipolar hard sphere fluid in C 3 . 
Here we will consider the 2D case which is a mere transposition so that we can skip many 
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details. Moreover we will consider only tinfoil boundary conditions, i.e. e = oo, which 
simplihes the algebra. The bare Green’s function is given by 

G„(r,. r,) ^ G„(r.,) = . (is) 

where '^(r) is the periodic Ewald potential. It satishes Poisson’s equation in C 2 

Aij{r) =-2n[Sc^{r) - , (19) 

where 


<^C2(r) 


=—^exp(ik-r), 

k 

is the periodical Dirac’s comb. Expanding "^(r) and Go(r) in Fourier series one hnds 

, , . 271 exp(ik ■ r) 

kT^O 


( 20 ) 


Go(r) = - ^ kk exp(ik • r) , 


(21a) 

(21b) 


kT^O 


where k = k/||k||. 

It shall proof useful in what follows to remark that one can rewrite the Ewald potential 
as 

71" o . 

( 22 ) 


■i/i(r) = - log r + , 

where 5D(r) is a harmonic function in the square and can thus be expressed quite generally 

as QQ 

OO 

S7p{r) = ^ cos{mip + a^) , (23) 

m=l 

where ip is the angle of r with the axis e^. and the constants a^, C(m are such that 'ip{r) is a 
periodical function. 

It is then easy to deduce from these prolegomena the two formulas 


TrGo(ri,r2) =j^ - 5c2(ri2) , 


(24a) 


111 

2ri2 ■ Go(ri, r2) ■ ri2 H- 5 “ + “ ~ l)amr™“^ cos(m(y9i2 + am) ■ (24b) 
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The computation of the dressed Green’s function from its definition (lllbp is conveniently 
made in Fourier space. Under the assumption of the locality of the dielectric tensor e(ri, r 2 ) 
one hnds the obvious result G = Go/e. Therefore Fulton’s relation flllap takes the explicit 
form 


X(ri, ra) = (e - l)I(ri, ra) + — -^Go(ri, ra). (25) 

e 

We stress that the above equation has been obtained under the assumption of the locality 
of the dielectric tensor e(ri,r 2 ). Therefore it should be valid only asymptotically, i.e. for 
points (ri, r 2 ) at a mutual distance ri 2 larger then the range ^ of e(ri, r 2 ). 

Taking the trace of Eq. fl25p . making use of Eq. fl24ap and integrating both ri and r 2 over 
the square C 2 one hnds the expression of the dielectric constant 


e - 1 = ^ < > , (26) 

where M = Mi is the total dipole moment of the square. 

We turn now our attention to the susceptibility tensor x(ri, r 2 ) which may be expressed 
in terms of the pair correlation function g{l, 2 ) where i = (rj,aj) {i = 1 , 2 ) denotes the 
position and the angle of dipole /Xj with axis e^,. One obtains that 

X(ri,r2) = |/I(ri,r2)+ 2j/py —j — h(l,2)siS2, (27) 

where y = 7r(3pfi‘^ and h = — 1 as usual. In the inhnite plane E 2 the pair correlation 

function g{l, 2 ) can be expanded on a complete set of rotational invariants among which the 
most important are 


- 00 ( 1 , 2 ) = 1 , 

(28a) 

^( 1 , 2 ) =Si ■ S 2 , 

(28b) 

-^(1, 2) =2(si ■ ri 2 )(s 2 ■ ri 2 ) — A(l, 2) , 

(28c) 


where ri 2 = ri 2 /ri 2 . 

In space C 2 the function g{l, 2) has the symmetry of the square and stricto sensu cannot 
be expanded onto these rotational invariants. However, following de Leeuw et a/. one 
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defines the projections 

M1.2)A(1,2) , (29a) 

h°(r,2) =2 ^ ^ h(l, 2)B(1,2) . (29b) 

(29c) 

Note that the two projections h^(ri2) and h^(ri2) are periodic functions which depend 
explicitely on the direction of vector ri2. The susceptibility tensor x(ri,r2) cannot be 
expressed in terms of these sole projections; however one can deduce from Eq. fl271) and the 
dehnitions (1291) the relations 


Trx(ri,r2) =2y +yph^iri 2 ) , (30a) 


2ri2 • x(ri,r2) ■ ri2 =yph^{ri 2 ) . 


(30b) 


The comparison of Eqs. (Elj), (IHOil . and (|25|1 yields the asymptotic behaviour of the projec¬ 
tions h^(ri2) and h^(ri2), he., for ||ri2|| > One has 


^asymp(^12) 


^fsymp(ri2) 
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yp^ in 2 J 

In actual simulations one rather computes angular averages of the functions h‘^{ri 2 ) and 
h^(ri2), he., 


A^™(r.2)= /'^A^™(ri2). 

Jo 271 

The asymptotic values of these averaged functions are simpler and given by 

(6-1)2 1 1 
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(33b) 


6 yp vrrja 

which are valid of course only for ^ < ri 2 < L/2. We note that, in the thermodynamic limit 
: he. for r hxed and L —>■ cx), one recovers the expected Euclidian behaviours {r) ~ 0 


(he. a short range function of r) and ~ (e — 1 

Euclidian plane E 2 without boundaries at inhnity (cf. Refs. 




Tiype) X l/r2 valid for the 

A- 
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Our last comment concerns Eq. fl2^ which can be recast as 


e -1 = 2/1 + 


d\ h^M 


( 34 ) 


f - 1 = i 1 + S / '>»(>■) Kt / '4,„p(n . 


(35) 


+2 

that we examine in the limit L ^ oo. We can then write 

'' ^ 4^rft^(r)|+^ [ 

E 2 J ^ JC 2 

where we have noted that, in the limit L —)■ oo, h‘^{r) —)■ h^{r) becomes an isotropic 
function. Making use of Eq. fl33al) to compute the second integral in one obtains 

(e-l )(6 + l) 


2 e 


= !/P + I 4^ h^(r) 


(36) 


This expression of e is precisely that obtained in space E 2 by various methods [3|, l4|. 


B. The sphere ^2 

We recall here the results of Ref. jd] for a fluid of bi-dipoles conhned on the surface of 
the sphere 1 S 2 . The dielectric constant is given by 
6-1 , ( 6 - 1)2 


26 


COS-00 = with 0 < ^jJo < 7r/2 , 


(37) 


where the fluctuation m2('?y)o) is given by 


m 2 (^y>o) 


TT/S/i^ 

R 


N N 

< 5 ^ 5 ^ Si ■ Sj Qd’o - Ai) > , 


(38) 


* J 

with S = 27 ri ?2 (surface of the northern hemisphere) and Q{x) the Heaviside step-function 
(0(x) = 0 for X < 0 and 0(x) = 1 for x > 0). In the MC simulations reported in this paper 
we retained the optimal choice 'ipo = tt/S. Asymptotically (he. for a large fixed r = R'lp S> ^ 
and '0 < ^/2), one has 


^fsymp(^) ~ 


(^- 1 )^ 1 1 
ype 271W' 1 + cosTp ’ 
(^- 1 )^ 1 1 
ype 2nR^ 1 — cos 27 


(39a) 

(39b) 


As for C 2 these asymptotic behaviours allow to recover from the formula fITTj) of the dielectric 
constant in space S 2 the expression fl36D in the thermodynamic limit. 
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IV. COMPARISONS OF THE TWO GEOMETRIES OF SIMULATION 


We performed standard MC simulations of the DHD fluid in the canonical ensemble 
with single particle displacement moves (translation and rotation) in both geometries 
and 1S2. Some elements of comparison are given in Table H] for three equilibrium states in 
the isotropic fluid phase of the model. We report values for the reduced internal energy 
per particle =< ^Udd > the contact values of the projections h^{a), and 

h^{a) of the pair correlation function g{l, 2 ), the compressibility factor Z = (iP/p ( P 
the pressure) with Z = Z^s + and Z^s = 1 + (7rp*/2)(7°°(cr), and the specihc heat 
Cy/kB = (< {(3UddY > — < PUdd >‘^)/N. As apparent in Table [H the agreement between 
the two methods of simulation is quite satisfactory. The values reported in the table were 
obtained for systems of ~ 1000 particles for which A"conf. = 5 — 10 x 10® conhgurations 
per particle were generated. Note that the hnite size scaling study of Ref. ^ gives, in the 
thermodynamic limit N ^ 00 , (3uoo = —1.79000(6) for state {p* = 0.7, p* = \/2) and 
(3uoo = —4.17138(15) for state {p* = 0.6, p* = 2), which shows that the data reported here 
are not very far from this limit. 

We have also tested the validity of the asymptotic behaviours of h'^(r), and h^{r) in both 
geometries. We display in Fig. [1] these functions as well as their asymptotic behaviours (l33|) 
and fl39l) for the state {p* = 0.6, p* = 2). The values of the dielectric constant which 
enter these asymptotic behaviours are those given in Table. [B As apparent on the hgures an 
excellent agreement between the MC data and the theoretical prediction is obtained. The 
small tails observed in h‘^{r) at large r, which differ signihcantly in the two geometries, are of 
primary importance to ensure that the dielectric constants e are identical in both geometries, 
within numerical uncertainties and hnite size effects, although given by completely different 
formulas. 


V. MC SIMULATIONS OF THE FLUID PHASE 

The homogeneous, isotropic huid phase is no more stable at low temperatures and com¬ 
plicated structures arise in this domain as indicated by some snapshots displayed in Fig. |2l 
At low densities, clusters of aligned dipoles, mostly organized into closed rings, appear at low 
temperatures and this topological structure becomes even more complex at higher densities. 
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In this low temperature regime the theory of the dielectric constant given in Sec. fllllj) be¬ 
comes incorrect and the predicted asymptotic behaviours of h^(r), and h^{r) are no more 
observed. Most probably the dielectric tensor, even if it exists, is no more isotropic and 
Fulton’s theory breaks down. In order to establish the thermodynamic stability of the high 
temperature phase we have followed the authors of Ref. jl4. 15| and computed the specific 
heat as a function of fi for some densities p* = 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, and 0.7. A 
peak in Ci,{p) should be a signal of the “transition” or the limit of stability of the fluid 
phase. Some curves Cy{p) are displayed in Fig. |3l They were obtained in the canonical 
ensemble for systems involving N ~ 1000 dipoles and runs of A^conf. = 5 — 10 x 10® con¬ 


figurations per particle. Table 


densities considered. In Ref 


ni provides the transition dipole moments for the different 
14| it has been pointed out that the polymerization transition 
may also be defined from the inflection point of ^ = Np/N as a function of dipole moment 
(or temperature T* = l//i*^) where Np is the number of particles belonging to a cluster. At 
the density p* = 0.05 where clusters are well defined we obtain a transition temperature in 
agreement with the value given in Table HIl 


VI. CONCLUSION 

In this paper we have studied the 2D DHD system by means of MC simulations performed 
either in a square with periodic boundary conditions or on the surface of a sphere. The 
interactions between dipoles have been chosen so as to satisfy the laws of electrostatics in 
the two geometries. With this precaution both methods lead to identical results for the 
thermodynamic, structural and dielectric properties of the system, at least for sufficiently 
large systems. A subtlety in the asymptotic behaviours of the pair correlation function, 
strongly depending on the geometry, has been predicted and observed in the MC experiments 
performed in the isotropic fluid phase. 

In the low temperature, low density part of the phase diagram a phase of living polymers 
of aligned dipoles organized into closed rings has been observed. At higher density the 
structure of this phase looks like an entangled structure of chains and rings. 

At these low temperatures the laws of macroscopic dielectrics seem to be violated. A 
polymerization transition line based on the maximum of the specific heat as a function of 
dipole moment is provided. The critical dipole moment p* at the transition from fluid to 
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polymeric phase increases slightly with density. 

We can contrast the present system with the one of 3D dipolar particles with centers 
of mass constrained to a monolayer or thin layer, at least if the dipoles are in-plane as 
it is the case at low temperatures. Such a quasi-two-dimensional (Q2D) system has been 
studied extensively in numerical simulations js, 16-19] in view of its relevance to various 


experimental situations. References to experimental works can be found in 


also Refs. 


[l5, 20 


16 


19| . see 


-l23|. 


In Q2D systems the head to tail interaction of two particles at contact is —2/i^ and 
antiparallel side by side interaction is —while in 2D the interactions of both types of 
arrangement are of similar strength —/i^. One would therefore expect that chaining is much 
favoured in the Q2D case. This is easily demonstrated by comparing structural properties 
obtained in simulations of both systems. Notwithstanding, the overall qualitative structural 
behaviour appears to be much the same at comparable densities (and short range interac¬ 
tion), especially at low temperature, i.e., formation of chains and rings. A notable difference 
between the Q2D and 2D systems is however that in the former system the spatial decay 
of the interaction is faster (l/r^) than the system dimension (D=2) i.e., of ’’short” range. 
Moreover, the angular dependence of the dipole-dipole interaction in Q2D systems is a linear 
combination of the 2D rotational invariants D{1,2) and A(l,2). 

Although 2D dipolar fluids do not exist per se in nature, the model could be used via 
various mappings for applications as, recently, for the hydrodynamics of two-dimensional 
microfluids of droplets. It is argued in Ref. [24 that droplet velocities show long-range 
orientational order decaying as 1/r^. 
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geometry p* 


I3u 

Z 

e 

0 

0 

h^(a) 

h^{a) 

Cv/ks 

C2 

0.7 

V 2 

-1.789 

3.990 

17.45 

4.346 

2.826 

4.578 

.712 

S2 

0.7 

V 2 

-1.790 

3.974 

: 17.83 

4.339 

2.814 

4.575 

.709 

C2 

0.6 

2.0 

-4.168 

1.680 

42.8 

5.144 

5.529 

7.131 

1.44 

S2 

0.6 

2.0 

-4.172 

1.662 

43.24 

: 5.132 

5.502 

7.117 

1.419 

C2 

0.1 

2.0 

-1.957 

0.53 

3.34 

9.48 

8.03 

15.91 

3.76 

S2 

0.1 

2.0 

-1.955 

0.53 

3.32 

9.45 

7.98 

15.85 

3.66 


TABLE I: Geometry of simulation, dimensionless numerical density p* , reduced dipole moment p* , 
reduced internal energy per particle I3u, dielectric constant e, compressibility factor Z = l3P/p, 
contact values h^{a), and h^{a) of some projections of the pair correlation function, and 

specific heat of the DHD fluid. For each state systems of = 1000 particles were considered in S 2 
and N = 1024 in C 2 . In both cases A^conf. = 5 — 10 x 10® configurations per particle were generated. 
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p* 0.05 0.1 0.2 0.3 0.4 0.5 0.6 

pI 2.35 ±0.15 2.4 ±0.1 2.5 ±0.1 2.5 ± 0.25 2.6 ±0.1 2.8 ±0.1 3.0 ±0.1 

TABLE II: pl{p*) for a system oi N = 1000 dipolar hard disks. 
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FIG. 1: Projections (Top ; C 2 , bottom ^ 2 ) and (Top ; ^ 2 , bottom C 2 ) for the state 
{p* = 0.6, p* = 2). Solid lines ; MC data, dashed lines : predicted asymptotic behaviours. In space 
C 2 ■ r < L/2 while in space S 2 ■ r < i?7r/2. 
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-20 -10 0 10 20 

FIG. 2: Snapshots of configurations of the DHD system in space C 2 . Top left : (p* = 0.1, = 2, 

N = 1024), Top right : {p* = 0.1, p* = 3, N = 1024), Bottom left : {p* = 0.6, p* = 2, N = 1024), 
Bottom right : {p* = 0.6, p* = 3.5, N = 2500). 
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